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Abstract 

We search for the lep ton-flavor- violating decay of the J/ip into an electron and a muon using (225.3 ± 
2.8) x 10 6 J/tp events collected with the BESIII detector at the BEPCII collider. Four candidate events 
are found in the signal region, consistent with background expectations. An upper limit on the branching 
fraction of S(J/V> -»• efi) < 1.5 x 1(T 7 (90% C.L.) is obtained. 

PACS numbers: 13.25.Gv, 11.30.Hv, 12.60.-i 
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I. INTRODUCTION 



With finite neutrino masses included, the Standard Model allows for Lepton Flavor Violation 
(LFV). Yet the smallness of these masses leads to a very large suppression, with predicted branch- 
ing fractions well beyond current experimental sensitivity. However, there are various theoretical 
models which may enhance LFV effects up to a detectable level. Examples of such model predic- 
tions, which often involve super- symmetry (SUSY), include SUSY-based grand unified theories 
(1, SUSY with a right-handed neutrino fl, gauge-mediated SUSY breaking 0|, SUSY with 
vector-like leptons [Hj], SUSY with R-parity violation [0], models with a Z' |@], and models violat- 
ing Lorentz invariance [@]. The detection of a LFV decay well above Standard Model expectations 
would be distinctive evidence for new physics. 

Experimentally, the search for LFV effects has been carried out using lepton (/x,r) decays, 
pseudoscalar meson (K,7r) decays, and vector meson (<j),J/ip,Y) decays, etc. For example, a recent 
search for the decay ofu + — > je + from the MEG Collaboration yields an upper limit of — > 
7e + ) < 2.4 x 1CT 12 [80, and in a similar search with r decays the BaBar Collaboration reports 
B(t + -> 7e+) < 3.3 x 10~ 8 §]. The latest results for neutral kaon and pion decays from the 
E871 Collaboration and the E865 Collaboration, respectively, are B{K\ — > fi + e~) < 4.7 x 10~ 12 
ifToh and B(n° fj, + e~) < 3.8 x 10" 10 iQ. The best cf> decay limit, based on 8.5 pb 1 of 
e + e~ annihilations at center-of-mass energies from y/s = 984 — 1060 MeV, is obtained by the 
SND Collaboration: B((f) — > (J> + e~) < 2.0 x 10~ 6 [Tj]. In the bottomonium system, based 
on about 20.8 million T(IS') events, 9.3 million T(25) events, and 5.9 million T(3S) events, the 
CLEOIII Collaboration presented the most stringent LFV upper limits, B(T(1S, 2S, 3S) — > fxr) < 
C(10~ 6 ) lfl3ll . For charmonium, the best limits come from the BESII Collaboration, who obtained 
BiJ/ip -> [ie) < 1.1 x 10" 6 Q, B(J/tp ->■ er) < 8.3 x 10" 6 and B(J/^ -> fir) < 2.0 x 10" 6 
[fTill from a sample of 58 million J/ip events. A recent sample of 225 million J/ip events [fl6ll 
collected with the much improved BESIII detector now allows for LFV searches in J '/ip decays 
with a significant improvement in sensitivity. We present here our results from a blind analysis of 
J/ip -> e ± /i T . 



II. BESIII DETECTOR AND MONTE CARLO SIMULATIONS 



The BESIII detector H17H at the BEPCII collider is a large solid-angle magnetic spectrometer 
with a geometrical acceptance of 93% of An solid angle consisting of four main components. The 
innermost is a small-cell, helium-based (40% He, 60% C3H8) main drift chamber (MDC) with 
43 layers providing an average single-hit resolution of 135 /im. The resulting charged-particle 
momentum resolution in our 1.0 T magnetic field is 0.5% at 1.0 GeV, and the resolution on the 
ionization energy loss information (dE/dx) is better than 6%. Next is a time-of-flight (TOF) sys- 
tem constructed of 5 cm thick plastic scintillators, with 176 detectors of 2.4 m length in two layers 
in the barrel and 96 fan-shaped detectors in the end-caps. The barrel (end-cap) time resolution of 
80 ps (110 ps) provides a 2a K/n separation for momenta up to 1.0 GeV. Continuing outward, 
we have an electromagnetic calorimeter (EMC) consisting of 6240 CsI(Tl) crystals in a cylindrical 
barrel structure and two end-caps. The energy resolution at 1.0 GeV is 2.5% (5%) and the position 
resolution is 6 mm (9 mm) in the barrel (end-caps). Finally, the muon counter (MUC) consisting 
of 1000 m 2 of Resistive Plate Chambers (RPCs) in nine barrel and eight end-cap layers; it provides 
a 2 cm position resolution. 

Our event selection and sensitivity, including backgrounds, are optimized through Monte Carlo 



(MC) simulation. The GEANT4-based simulation software BOOST [18] incorporates the geom- 
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etry implementation simulations and material composition of the BESIII detector, the detector 
response and digitization models as well as the tracking of the detector running conditions and 
performances. The generic simulated events are generated by e + e~ annihilation into a J/tp meson 
using the generator KKMC lfl9ll at energies around the center-of-mass energy y^s = 3.097 GeV. 
The beam energy and its energy spread are set according to measurements of BEPCII, and initial 
state radiation (ISR) is implemented in the J/ij) generation. The decays of the J/ip resonance 
are generated by EVTGEN ||20|l for the known modes with branching fractions according to the 
world-average values l|21|l . and by lundcharm Il22ll for the remaining unknown decay modes. 



HI. EVENT SELECTION 

We search for events in which J/ ip decays into an electron and a muon. Candidate signal events 
are required to have two well-measured tracks with | cos#| < 0.8 and zero net charge, consistent 
with originating from the interaction point. Here, 9 is the polar angle with respect to the beam axis 
and the closest approach of each track to the interaction point must be less then 5 cm (1 cm) in the 
beam direction (in the plane perpendicular to the beam). To reject cosmic rays, the TOF difference 
between the charged tracks must be less than 1.0 ns. The acollinearity and acoplanarity angles 
between two charged tracks are required to be less than 0.9° and 1.4°, respectively, to reduce other 
backgrounds. 

In order to suppress the radiative events from e + e~ — > je + e~ and e + e~ — > 7/i + /i~, we veto 
events with one or more good photon candidates passing the following requirements. Candi- 
date showers reconstructed in both the EMC barrel region (| cos#| < 0.8) and in the end-caps 
(0.86 < | cos 6*1 < 0.92) must have a minimum energy of 15 MeV. Showers in the angular range 
between the barrel and end-cap are poorly reconstructed and are not considered. Showers caused 
by charged particles are eliminated by requiring candidates to be more than 20 degrees away from 
the extrapolated positions of all charged tracks. Requirements on EMC cluster timing suppress 
both electronic noise and energy deposits unrelated to the event. 

The above selection criteria retain events with back-to-back charged tracks and no obvious extra 
EMC activity. Most of the remaining events originate from the background processes J/ip — > 
e + e~,J/vfj — > J /if) — > ir + n~ , J/vp — > K + K~,e + e~ — > e + e~( / y) and e + e~ — > /i + /i"(7), 

In order to suppress the these background events, we identify electrons and muons based on the 
information of the MDC, EMC and MUC sub-detectors. The requirements are determined using 
electron, muon, pion and kaon samples from J/ip — > e + e~ , ,7r + n~ , K + K~ MC events. 
Electron identification requires no associated hits in the MUC and 0.95 < E/p < 1.50, where 
E is the energy deposited in the EMC and p is the momentum measured by the MDC. Also, 
the absolute value of XdE/dx fr° m comparing the dE/dx measurement with the expected electron 
signal should be less than 1.8. Fig. Q] shows the E/p and XdE/dx distributions for electrons, which 
are well-separated from other particles. Muon identification uses the barrel MUC system which 
covers | cos^l < 0.75. Charged tracks are required to have E/p < 0.5 and a deposited energy in 
the EMC 0.1 < E < 0.3 GeV. We require the penetration depth in the MUC to be larger than 40 
cm; if the track penetrates more than three detecting layers in the MUC, we also require the MUC 
track fit to have x 2 < 100. Finally, the XdE/dx va l ue from the dE/dx measurement calculated 
with the electron hypothesis must be less than —1.8. The simulated distributions of the deposited 
energy in the EMC and the penetration depth in the MUC are shown in Fig. [2l 

Our final analysis of event yields for J/ip — > e + [T~ is performed with the two variables | T,p\/y/s 
and E vis / ^fs, where |Ep| is the magnitude of the vector sum of the momentum in the event, E vis 
is the total reconstructed energy (calculated using \/p 2 + m 2 with each track momentum p), and 
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FIG. 1: The distributions of (a) E/p and (b) XdE/dx ^ or tne simulated electron, muon, pion and kaon 
samples. 
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FIG. 2: The distributions of (a) the penetration depth in the MUC and (b) the deposited energy in the EMC 
for the simulated muon, pion and kaon samples. 

y/s is the center-of-mass energy. Our signal region is defined by 0.93 < E vis /y/s < 1.10 and 
\^p\/y/s ^ 0-10' which correspond in each case to about two standard deviations as determined 
by MC simulation. 

The analysis is done in a blind fashion in order not to bias our choice of selection criteria. Be- 
fore examining the signal region, all selection criteria were optimized based on simulated samples 



6 



with a sensitivity figure-of-merit (FOM) defined as the average upper limit from an ensemble of 
experiments with the expected background and no signal, 



FOM 



T,n„ =0 P{Nobs\N exp ) ■ UL{N obs \N t 



exp ) 



(1) 



where e is the detection efficiency determined with a sample of 100,000 simulated J/ip — > e/i 
events, N exp is the expected number of background events based on background process simula- 
tions, N b s is the number of observed candidate events, P is the Poisson probability, and UL is the 
upper limit on the signal calculated with the Feldman-Cousins method at 90% C.L. [230. hi addi- 
tion to the signal MC samples, six background MC samples, each with twice the statistics of the 
data sample, are employed to optimize the selection criteria: J/ip — > e + e~, J/ip — > J/ip — > 

7r + 7r - , J/ip — > K + K~, e + e~ — > e + e~(^), and e + e~ — > /i + /i~(7). 

After applying the optimized selections criteria, four candidate events remain in our signal 
region, see Fig. [3] The detection efficiency for signal is determined to be (18.99 ± 0.12)%. Using 
an inclusive sample of simulated J / ip decays with four times the size of our data sample, we find 
nineteen background events surviving in the signal region. This yields a predicted background of 
^ = (4.75=1:1.09). 




0.9 1 1.1 

E v i s / \fs 
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FIG. 3: A scatter plot of E v i s / yfs versus jSpj/y^ for the J/ip data. The indicated signal region is defined 
as 0.93 < E vis /s/s < 1.10 and \Ep\/y/s < 0.1. 



IV. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties originate from imperfect knowledge of efficiencies for the electron 
and muon tracking requirements electron and muon identification, acollinearity and acoplanarity 
requirements, the photon candidate veto, and the number of J/ip events. 
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A. Tracking efficiency 



Control samples of ijj' — > n + n^J/ip, J/ijj — >■ e + e~, selected from 106 M ijj' data events 

and 106 M ip' inclusive MC events are used to study the possible differences in the tracking effi- 
ciency between data and MC events. To determine the tracking efficiency of electrons, we select 
events with at least three charged tracks. Two tracks with low momentum, p < 0.5 GeV/c, and 
with opposite charge are interpreted as the pions. After requiring the recoiling mass opposite these 
two pions to satisfy \M re f™± - 3.097| < 10 MeV, we obtain ip' -> 7T + 7r J/ip candidates. For the 
e~ selection, at least one track is required to have a negative charge, a momentum in the region 
from 1 .0 GeV/c to 2.0 GeV/c, and a deposited energy in the EMC greater than 1 .0 GeV. With these 
three tagging tracks, n + , ix~ , and e~, the total number of e + tracks, N°+, can be determined by 
fitting the distribution of mass recoiling from the 7r + 7r~e~ system, M^ e( ^ il e . In addition, one can 
obtain the number of detected e + tracks, N*+, by fitting M^ ec ^ u e , after requiring all four charged 
tracks to be reconstructed. The tracking efficiency of e + is then obtained as e e + = N^ + /N® + . 

Similarly, we can obtain the tracking efficiency for e~, /i + and fi~. The difference between 
data and MC simulation is found to be about 1 .0% in each of these four cases, which is taken as a 
systematic uncertainty for tracking. 

B. Particle identification 

Clean samples of J/ip(e + e~) — > e + e~ with backgrounds less than 1% selected from data and 
inclusive MC events are employed to estimate the uncertainty of the e ± identification. The event 
selection criteria for this control sample are identical to those for our signal channel, including 
two good charged tracks and no good photon. The track with higher momentum is required to 
satisfy the e ± identification criteria described previously, and the other track is used for the e T 
identification study. 

The electron identification efficiency, obtained by comparing the number of events with and 
without electron identification criteria applied on the selected control sample, is defined by: 
e PiD = N evt (w/ PID)/N evt (w/o PID), where N evt is the number of events extracted from 
the control sample. It is found that the average efficiency difference between data and MC simu- 
lation is 0.62% for the track momentum range 1.4 — 1.7 GeV/c, which is taken as the systematic 
error for electron identification. Applying a similar method, we study the systematic error of the 
^ identification using the control sample J /ip(e + e~) — > /i + yU~. We apply corrections based on 
data-MC differences, and a residual uncertainty of 0.04% is obtained for the muon identification 
in the momentum range 1.4 — 1.7 GeV/c. 

C. Acollinearity and acoplanarity angles 

A control sample of J/ij) — > is employed to estimate the uncertainty due to the acollinear- 

ity and acoplanarity angle requirements. We obtain the corresponding selection efficiency by com- 
paring the number of events with and without imposing the acollinearity and acoplanarity angle 
requirements on the the selected control sample. We find an efficiency difference between data 
and MC simulation of 2.83%, which is taken as a systematic uncertainty for acollinearity and 
acoplanarity angle requirements. 
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D. Photon veto 



We expect no good photon candidates to be present in J/ijj — > , and therefore choose this 
channel as a suitable control sample. The event selection criteria for this control sample are similar 
to those in sub-section B. By comparing the numbers of events before and after imposing the 7- 
veto criteria on the selected control sample, we can obtain the corresponding selection efficiency. 
We find that the difference in efficiency between data and MC simulation is 1 .19%, which is taken 
as a systematic uncertainty for photon veto. 



The uncertainty in the number of J/ij) is 1.24% [16]. Table U summarizes the systematic er- 
ror contributions from different sources and the total systematic error is the sum of individual 
contributions added in quadrature. 



TABLE I: Summary of systematic uncertainties (%). 



Sources Error 
tracking 1 .00 

^ tracking 1.00 
e ± ID 0.62 
/i ± ID 0.04 
Acollinearity, acoplanarity 2.83 
Photon veto 1.19 

Nj^ k24_ 

Total 3.65 



V. RESULTS 



We observe four candidate events with an expected background of 4.75 ± 1.09, and therefore 
set an upper limit on the branching fraction of J/tp — > e/i, based on the Feldman-Cousins method 
with systematic uncertainties included. The upper limit on the number of observed signal events 



at 90% C.L., N%£, of 6.15 if obtained with the POLE program Il25ll . Here, the number of expected 
background events, the number of observed events, and the systematic uncertainty are used as the 
input parameters. The upper limit on the branching fraction is given by 

N U , L 

B{J/^ e/i) < —2*2-, (2) 
J/,;, ■ e 

where Njm is the total number of J/ip events, and e is the detection efficiency. Combining, we 
find a 90% C.L. upper limit on the branching fraction of B(J/ip — > e/i) < 1.5 x 10~ 7 . 



VI. SUMMARY 



Using 225.3 ± 2.8 x 10 6 J/tp events collected with the BESIII detector, we have performed 
a blind analysis searching for the lepton flavor violation process J/ip — > efi. We observe four 
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candidate events, consistent with our background expectation. The resulting upper limit on the 
branching fraction, B{J/ip — > e/x) < 1.5 x 10~ 7 (90% C.L.), is the most stringent limit obtained 
thus far for a LFV effect in the heavy quarkonium system. 
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